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Present-day implementations of quantum information processing rely on two 
widely different types of quantum bits (qubits). On the one hand, microscopic 
systems such as atoms or spins are naturally well decoupled from their environ- 
ment and as such can reach extremely long coherence times [H [2]; on the other 
hand, more macroscopic objects such as superconducting circuits are strongly 
coupled to electromagnetic fields, making them easy to entangle [31 H] although 
with shorter coherence times pl|6j. It thus seems appealing to combine the two 
types of systems in hybrid structures that could possibly take the best of both 
worlds. Here we report the first experimental realization of a hybrid quantum 
circuit in which a superconducting qubit of the transmon type [H [7] is coher- 
ently coupled to a spin ensemble consisting of nitrogen-vacancy (NV) centers 
in a diamond crystal ^Sj via a frequency-tunable superconducting resonator [9] 
acting as a quantum bus. Using this circuit, we prepare arbitrary superpositions 
of the qubit states that we store into collective excitations of the spin ensemble 
and retrieve back later on into the qubit. We demonstrate that this process 
preserves quantum coherence by performing quantum state tomography of the 
qubit. These results constitute a first proof of concept of spin-ensemble based 
quantum memory for superconducting qubits |10H12| . As a landmark of the 
successful marriage between a superconducting qubit and electronic spins, we 
detect with the qubit the hyperfine structure of the NV center. 

Superconducting qubits have been successfully coupled to electromagnetic [13] as well as 
mechanical [14j resonators; but coupling them to microscopic systems in a controlled way 
has up to now remained an elusive perspective - even though qubits sometimes turn out to be 
coupled to unknown and uncontrolled microscopic degrees of freedom with relatively short 
coherence times [15j. Whereas the coupling constant g of one individual microscopic system 
to a superconducting circuit is usually too weak for quantum information applications, en- 
sembles of such systems are coupled with a constant g^/N enhanced by collective effects. 
This makes possible to reach a regime of strong coupling between one collective variable 
of the ensemble and the circuit. This collective variable, which behaves in the low excita- 
tion limit as a harmonic oscillator, has been proposed [ 10l - [T2] as a quantum memory for 
storing the state of superconducting qubits. Experimentally, the strong coupling between 
an ensemble of electronic spins and a superconducting resonator has been demonstrated 
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spectroscopically p!6l - [T8] . and the storage of a microwave field into collective excitations of 
a spin ensemble has been observed very recently [T9l [20] . These experiments were however 
carried out in a classical regime since the resonator and spin ensemble behaved as two cou- 
pled harmonic oscillators driven by large microwave fields. In the perspective of building 
a quantum memory, it is instead necessary to perform experiments at the level of a single 
quantum of excitation. For that purpose, we integrate for the first time on the same chip 
three different quantum systems : an ensemble of electronic spins, a superconducting qubit, 
and a resonator acting as a quantum bus between the qubit and the spins. A sketch of the 
experiment is shown in Fig. [TJ 

The spin ensemble NV consists of ^ lO"*^^ negatively-charged NV color centers [8j in a 
diamond crystal. These centers have an electronic spin S = 1^ with electron spin resonance 
(ESR) transition frequencies (jU± between energy levels ms = and ms = ±1 of about 
2.88 GHz in zero magnetic field (see Fig. The electronic spin of the NV center is 

further coupled by hyperfine (HF) interaction to the spin-one ^^A^ nucleus, which splits (jJ± 
into three peaks separated by 2.2 MHz [2T] (see Supplementary Material). In our experiment, 
the diamond crystal is glued on top of the chip, and the degeneracy between states ms = ±1 
is lifted with a B^y = 1.4 mT magnetic field applied parallel to the chip and along the [1,1,1] 
crystalline axis. The NV frequencies being sensitive only to the projection of Bjyy along the 
N — V axis, two groups of NVs thus experience different Zeeman effects: those along [1,1,1] 
(denoted /) and those along either of the three other (1,1,1) axes (denoted /// as they are 
3 times more numerous). This results in four different ESR frequencies uj±i^±in. 

The qubit Q is a Cooper-pair box of the transmon type [5l [7] with transition frequency ujq 
between its ground state l^') and excited state |e). It is coupled to a nonlinear resonator R 
which is used to read-out its state, as in related circuit quantum electrodynamics experiments 
p2j. Single-qubit rotations are realized by applying microwave pulses at ujq through R. 
Qubit state readout is performed by measuring the phase of a microwave pulse refiected on 
i?, which depends on the qubit state; the probability Pg to find the qubit in its excited state 
is then determined by repeating ^ 10^ times the same experimental sequence. 

The quantum bus S, a superconducting resonator with quality factor ^ 10^, is electro- 
statically coupled to the qubit and magnetically coupled to the spin ensemble. In order to 
bridge the difference in frequency between Q and NV^ the bus frequency ujb can be tuned 
on a nanosecond time scale [23] by applying current pulses through an on-chip line, inducing 
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a magnetic flux $ through a SQUID embedded in S [9]. More information on the qubit 
readout and setup can be found in the Supplementary Material. 

We flrst characterize our hybrid circuit by spectroscopic measurements. The NV fre- 
quencies and coupling constants are obtained by measuring the microwave transmission 
|5'2i(cl;)| through the bus, while scanning its frequency (jOb{^) across the NV resonance. 
Vacuum Rabi splittings are observed when ujb matches the spin resonance frequency at 
cj+//27r = 2.91 GHz, cj_//27r = 2.84 GHz, cj+////27r = 2.89 GHz, and cj-////27r = 2.865 GHz 
(see Fig. [T]). From the data we extract the coupling constants g±i/27r = 2.9 MHz and 
g±iii/'27T = 3.8 MHz, the difference between the two values resulting essentially from the 
larger number of NV centers in group ///. Qubit spectroscopy is performed by scanning 
the frequency of a microwave pulse applied through i?, and by measuring Pg, which yields 
(jOq /27r = 2.607 GHz. This spectroscopy, measured while scanning ujb across cjg, shows an 
anticrossing (see Fig. ^) that yields the coupling constant gq/^n = 7.2 MHz between Q 
and B. 

Throughout the experiments reported in the following, the spins and qubit frequencies are 
kept flxed, and only ujb is varied in order to transfer coherently quantum information between 
Q and NV. For this purpose, a key operation is the qubit-bus SWAP gate that transfers 
an arbitrary qubit state a\g) + /3 \e) into the corresponding photonic state a\0)^ + /3\1)^ 
of the bus, leaving the qubit in l^'). This SWAP gate could be performed by tuning ujb 
in resonance with ujq for a duration 7T/2gQ [24]. Here we prefer instead to adiabatically 
sweep ujb across ujq as this sequence is more immune to flux noise in the SQUID loop of B 
(see Supplementary Material). This adiabatic SWAP (aSWAP) achieves the same quantum 
operation as the resonant SWAP except for an irrelevant dynamical phase. The experiments 
then proceed by combining single-qubit rotations, aSWAP gates, and flux pulses placing B 
and NV in and out of resonance for properly chosen interaction times. 

We apply such a sequence with the qubit initially prepared in |e) (see Fig. [2]). A flrst 
aSWAP converts |e) into the bus Fock state B is brought in or near resonance with the 
spin ensemble for a duration r; the resulting B state is then transferred back into the qubit, 
which is flnally read-out. Figure [2j3 shows the resulting curve Pei^) when the bus is brought 
in resonance either with uj-ni or uj^j. An oscillation in P^ is observed, revealing a storage in 
the spin ensemble of the single quantum of excitation initially in the qubit at Tsjh = 65 ns 
or ^s,i = 97 ns, and a retrieval back into the qubit at Trjn = 116 ns or Trj = 146 ns. The 
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Figure 1: Description of the hybrid quantum circuit demonstrated in this work. a,b Three- 
dimensional sketch of the device and corresponding electrical scheme. The ensemble NV of elec- 
tronic spins (magenta) consists of IQ-*^^ NV centers in a diamond crystal glued on the chip surface. 
The transmon qubit Q (in red) is capacitively coupled to a resonator R (in blue) made nonlinear 
with a Josephson junction and used to read-out its state. The bus B (in yellow) is electrostatically 
coupled to Q and magnetically coupled to NV. B contains a SQUID that makes its frequency 
(jJb{^) tunable by changing the flux $ in the SQUID loop applied via a fast on-chip current line 
(in green). A magnetic field B]\iy is applied parallel to the [1, 1, 1] crystallographic axis, c, (lower 
left inset) Energy level structure of NV centers. Transitions between ms = and ms = ±1 at 
frequency (jJ± are further split in three resonance lines due to the hyperfine interaction with the ^^N 
nuclear spin [21j. (main panel) Two-dimensional plot of the transmission 15211 (^5^) through B in 
dB units, with $ expressed in units of the superconducting flux quantum $0 = /i/2e, for a fleld 
Bnv = 1.4 mT applied to the spins. Color scale goes from —55 dB (green) to —30 dB (magenta). 
Four vacuum Rabi splittings are observed whenever cub matches one NV center resonance frequency. 
The four frequencies correspond to the uj±ijn transition frequencies of one of two distinct families 
of NV centers, being either along the [1, 1, 1] crystal direction parallel to Bjsiv {!•, in red), or along 
one of the three other possible (1,1,1) axes (labelled ///, in blue), as shown in upper right panel, 
(main panel, bottom right) Qubit excited state probability Pg as a function of the frequency of the 
exciting microwave and Color scale goes from 0.1 (purple) to 0.3 (yellow). When ujb matches 
the qubit frequency ujq — 2.607 GHz, the qubit spectrum shows an anticrossing demonstrating its 
coupling to B with constant ^g/27r = 7.2 MHz. 
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Interaction time,T (ns) Interaction time.x (ns) 

Figure 2: Storage and retrieval of a single quantum of excitation from the qubit to the spin ensemble, 
a, Experimental sequence showing the microwave pulses used for exciting the qubit in |e) (red) and 
for reading it out (blue), as well as transition frequencies of the quantum bus (orange), qubit (red), 
and spins (magenta), b. Experimental (red dots) and theoretical (black line - see text) probability 
Pe{T) for CjJb tuned to uj-m (top graph) or uj-i (bottom graph), showing the storage and retrieval 
times Ts and r^. c. Two-dimensional plot of Pg versus interaction time r and flux pulse height 
showing resonance with the four spin groups. Chevron-like patterns are observed, showing a faster 
oscillation with reduced amplitude when cj^ is detuned from the spin resonance, as expected. Note 
that the difference between the uj- and ujj^ patterns in the same NV group is simply caused by the 
non-linear dependence of cj^ on $ 

fidelity of this storage-retrieval process, defined as Pe(r^)/Pe(0), is 0.14 for group /// and 
0.07 for group /. These relatively low values are not due to a short spin dephasing time, but 
rather to an interference effect caused by the HF structure of NV centers, as evidenced by the 
non-exponential damping observed in Pei'r). These measurements are accurately reproduced 
by a full calculation of the spin-resonator dynamics (see Suppl. Mat. and [2211251126]) taking 
into account this HF structure, with the line width of each HF peak as the only adjustable 
parameter. A linewidth of 1.6 MHz is in this way determined for the spins in group /, and of 
2.4 MHz for group ///, this larger value being likely due to a residual misalignment of Bjsiv 
from the [1,1,1] crystalline axis causing each of the three < 1,1,1 > N — V orientations 
non-coUinear with the field to experience slightly different Zeeman shifts. We finally note 
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Figure 3: Storage and retrieval of a coherent superposition of states from the qubit to the spin 
ensemble, a, Experimental pulse sequence: the qubit is prepared by a 7r/2 pulse in state (|^) + 
|e))/\/2 , which is transferred to B by an aSWAP. B is then immediately tuned to a;_i/27r = 
2.84 GHz for an interaction time r. The quantum state of B is then transferred back to the qubit 
by an aSWAP. Quantum state tomography is finally performed to determine the qubit state by 
applying either /, X, or Y operation to the qubit. b, Trajectory of the qubit Bloch vector on the 
Bloch sphere (bottom inset), and its projection on the equatorial plane, c, Modulus and phase of 
the off-diagonal element pge of the qubit density matrix as a function of interaction time r. 

that in both curves shown in Fig. ^ Pe{^) tends towards 0.08 at long times, as is also found 
with the qubit initially in l^'). This proves that the collective spin variable coupled to B is, 
as requested for experiments in the quantum regime, in its ground state |0)_j_jjj with a 
large probability ^ 0.92 at equilibrium, which corresponds to a temperature of ^ 50 mK. 
Varying both ujb and r with the same pulse sequence, we observe similar storage-retrieval 
cycles at all four spin frequencies (see Fig. 

In addition to storing a single excitation from the qubit, one has to test if a coherent 
superposition of states can be transferred to the spin ensemble and retrieved. For that, we 
now perform the aSWAP and bring ujb in resonance with uj^i after having initialized the 
qubit in d^') + \e))/^/2 instead of |e), and we reconstruct the Bloch vector of the qubit by 
quantum state tomography at the end of the sequence. More precisely, we measure (ctx), 
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Figure 4: Ramsey-like experiment on the spin ensemble at the single-photon level, a, Experimental 
pulse sequence: the qubit is prepared in its excited state |e) by a tt pulse; the state |e, 0^) is then 
adiabatically transferred to \g, 1b) by an aSWAP. A fast flux pulse subsequently brings cub onto 
cj_i, and then lets B and the spins from group —I interact for half a swap time Tsj/2, generating an 
entangled state of the two systems. B is then detuned from the spins by Suj /27r = 38 MHz during 
a time r, and a second half-swap is performed. The quantum state of B is then transferred back to 
the qubit, which is finally read-out. b, Measured (red circles) and calculated (black line - see text) 
probability Pe{^)j as well as its Fourier transform (inset) revealing the NV centers HF structure. 

(cTy) and {az) by using 7r/2 rotations around y, X, or no rotation at all (/) prior to qubit 
readout. After substracting a trivial rotation around Z occurring at frequency {uj^j — ujq)^ 
we reconstruct the trajectory of this Bloch vector as a function of the interaction time r. It is 
plotted in Fig. |3| together with the oflF-diagonal element Pg^ of the final qubit density matrix, 
which quantifies its coherence. We find that no coherence is left in the qubit at the end of 
the sequence for r = r^j, as expected for a full storage of the initial state into the ensemble. 
Then, coherence is retrieved at r = r^j^ although with an amplitude ^ 5 times smaller 
than its value at r = (i.e. without interaction with the spins). Note the tt phase shift 
occurring after each storage-retrieval cycle, characteristic of 27r rotations in the two-level 
space {|1^,0_/) , |0b, I-/)}- The combination of the results of Figs. [2] and |3] demonstrates 
that arbitrary superpositions of the two qubit states can be stored and retrieved in a spin 
ensemble - although with limited fidelity - and thus represents a first proof-of-concept of a 
spin-based quantum memory for superconducting qubits. 

To evaluate the time during which quantum coherence can be stored in the ensemble. 
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we perform a Ramsey-like experiment on the spin ensemble at the single-photon level (see 
Fig. [i]): we initially prepare the qubit in |e), transfer its state to 5, then tune ujb to uj^j 
for a duration = rg^_//2, after which ujb is suddenly detuned by 6(jj/27t = 38 MHz 
for a time r. At this point, the joint bus-spin ensemble state is an entangled state 
(|1^, 0_/) + e^"^ |0b, 1-/)) /a/2 with a phase ip = 5ujt. B is then put back in resonance with 
the spins for a second interaction of duration Tt^i2 that converts the phase into population 
of |1^, 0_/). This population is finally transferred to the qubit, and read-out. Oscillations at 
frequency 5uj are observed in Pe{T) as seen in Fig. |4| confirming that the resonator and the 
spins are entangled after the first 7r/2 pulse. These oscillations are modulated by a beating 
pattern, with an overall damping of the oscillations envelope in ^ 200 ns. Quite remark- 
ably, this beating observed in the qubit excited state probability is directly caused by the 
HF structure of NV centers, as proved by the Fourier transform of Pe{T) which shows the 
three HF lines. The full calculation of the system dynamics quantitatively captures both the 
beatings and the oscillations damping, which is thus completely explained by the 1.6 MHz 
inhomogeneous linewidth of each HF line taken into account in the theory. 

The previous results suggest that the storage of quantum information in the NV cen- 
ters ensemble is at present limited both by its HF structure and by the inhomogeneous 
broadening of its resonance. This broadening is attributed to dipolar interactions between 
the NV centers and residual paramagnetic impurities (likely neutral nitrogen atoms) in the 
diamond crystal. Purer crystals could thus greatly improve the present performance of our 
device. Note that the hyperfine coupling to the nuclear spin of ^^A^ could be turned into a 
useful resource if quantum information was transferred from the electron spin to the nuclear 
spin degree of freedom, which has much narrower linewidth. Finally refocusing techniques 
borrowed from quantum memories in the optical domain [27j should also lead to increase in 
the storage time by two orders of magnitude. 

In conclusion our experiments bring the first proof of concept of a spin-based quantum 
memory for superconducting qubits. In a longer-term perspective, they open the way to 
the implementation of genuine quantum lab-on-chips, where superconducting qubits would 
coherently interact with electron and nuclear spins as well as optical photons. 

Note: During redaction of this manuscript we became aware of related work demonstrat- 
ing the coherent dynamics of a fiux qubit coupled to an ensemble of NV centers in diamond 
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Methods Summary 

Superconducting circuit parameters The transmon parameters are measured by 
standard spectroscopy, yielding a Josephson energy Ej/h = 5.2 GHz and a Coulomb energy 
for a Cooper-pair Ec/h = 0.66 GHz. Its relaxation time Ti = 1.75 /xs and coherence time 
T2 = 2.2 /is at the bias point used in this work were measured by standard pulse sequences. 

The bus resonator B could be tuned from a maximum frequency cc;^(0)/27r = 3.004 GHz 
down to 2.5 GHz. Its quality factor is 2 • 10^ at $ = in presence of the diamond sample 
(which thus does not introduce dielectric losses contrary to what was reported in p!6l [18] ). 
This quality factor degrades progressively as cj^ is tuned towards lower frequencies. However 
the bus resonator energy relaxation time Tcav = 1-5 /xs was measured using the qubit as 
explained in [29j, and was found not to depend on the flux bias. This indicates that the 
quality factor degradation is due to low-frequency noise, likely flux noise in the SQUID loop. 

The readout resonator R has a frequency ujr/2ti = 3.468 GHz and a quality factor Q = 
500. Its nonlinearity is brought by a Josephson junction of critical current 650 nA yielding 
a Kerr constant K/ujr = —4.5 • 10~^ [30]. Readout pulses have a frequency 3.456 GHz. The 
coupling between Q and R is estimated to be gR/27T = 30 MHz. 

Theory Each spin is modelled by an effective harmonic oscillator bj of frequency cjj, 
following the Holstein-Primakoff approximation valid in the low-excitation limit. The qubit, 
bus resonator, and spin ensemble are described by Hamiltonians Hq/K = —{ujq/2)(Jz-> 
Hb/^i = ujB{^)ci^ci and Hnv/J^ = being the Pauli matrix, a being the bus 

resonator annihilation operator. Coupling of the qubit to the bus resonator is described by 
a Jaynes-Cummings Hamiltonian Hq^b/^ = Qqi^^^ci + h.c) where a+ is the qubit raising 
operator. Coupling between the bus resonator and the spin ensemble is described by a Tavis- 
Cummings Hamiltonian Hb-nv = X^jLi gj{b^jCi + h.c.). This Hamiltonian can be rewritten 
as Hb-nv/^ = '}2k=-iii -I +1 +111 9k (^^^k^ + h.c^ where gx represents the collective cou- 
pling between B and each of the four spin groups {±I^±III) and bx = {^/ 9k) "^2^=1 9 i^i 
is the collective excitation annihilation operator. Excited states of the spin ensemble are 
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defined for each family by applying collective operators bx and 6^ to the ground state \0k) 
(for instance = b^_j |0_/)). The spin-resonator dynamics can be calculated from this 

model, as explained in [25l [26] and in the Supplementary Material. 

Adiabatic pulse parameters Our adiabatic SWAP operation proceeds as follows: ujb 
starts at 2.52 GHz, is first ramped up to 2.589 GHz in 60 ns, then to 2.643 GHz in 350 ns, 
then to 2.687 GHz in 40 ns. See the Supplementary Information for more details on the 
pulse optimization. 
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Supplementary information for "Hybrid quantum circuit with a superconducting 
qubit coupled to a spin ensemble" 

A. Diamond sample preparation and characterization 

The sample we use is a poHshed (110) plate of dimensions 2.2 x 1 x 0.5 mm^ taken from 
a type-lb HPHT crystal which contained 40 ppm of neutral substitutional nitrogen (the PI 
centers) as measured by IR absorption. It has been irradiated by 2 MeV electrons at 700° C 
with a total dose of 6.4 x lO^^e/cm^ and annealed at 1000° C for 2 hours in vacuum. The 
high temperature irradiation was employed to minimize the residual unwanted defects. 

The resulting concentration of negatively-charged NV centers was measured by comparing 
the sample photoluminescence (PL) to the one obtained from a single NV center. For that 
purpose, a continuous laser source operating at 532 nm wavelength was tightly focused on 
the sample through a high numerical aperture oil-immersion microscope objective. The NV 
center PL was collected by the same objective, spectrally filtered from the remaining pump 
light and directed to a silicon avalanche photodiode. After calibration of the PL response 
associated with a single NV center, a PL raster scan of the sample directly indicates the NV 
center content since the excitation volume is known. As shown in Fig. [5](a), the NV center 
concentration is found rather inhomogeneous over the sample, with an average density of 
4.4 X 10^'' cm~^. Such a value is in good agreement with the spin ensemble-resonator coupling 
constants reported in the main text of the manuscript. 
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Figure 5: Characterization of the diamond sample at room temperature, a, Two- 
dimensional map of the NV center density for a depth of 5 jj^m inside the diamond sample, b, 
Optically detected ESR spectrum with a static magnetic field B ^ 1.5 mT applied along the [111] 
axis of the diamond crystal. Four ESR lines are observed, corresponding to rris = ^ = ±1 
spin transitions for the two subsets of NV centers crystalline orientations (±111 and ±1). c, Hy- 
perfine structure. Data fitting with Lorentzian functions leads to Ahf = 2.17 =b 0.02 MHz as 
expected from hyperfine interaction with ^^N nuclear spins, d, Ramsey fringes recorded for the 
subset of [lll]-oriented NV centers with a microwave detuning 5 = 13 MHz from the ESR line at 
(jj^j /27r = 2.915 GHz. The laser pulses used in the Ramsey sequence have a duration of 30 fis and the 
laser power is set at 40 mW. For data analysis, the NV center PL recorded during the first 10 /xs of the 
laser pulses is used for spin-state read-out while the PL recorded during the last 10 fis is used as ref- 
erence. The red solid line is data fitting with the function exp[— r/TI] x Ya=-i [27r(5 + iAhf)T]. 
The inset shows the Fourier-transform of the free induction decay. Solid lines are data fitting with 
Lorentzian functions, e. Measurement of the coherence time T2 for the subset of [lll]-oriented NV 
centers using a7r/2 — r — tt — r — 7r/2 spin echo sequence. Data fitting with an exponential decay 
leads to T2 = 7.3 ± 0.4 fis. 
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Electron spin resonance (ESR) measurements were performed at room temperature by 
applying a microwave field through a copper microwire directly spanned on the diamond 
surface. In addition, a static magnetic field B ^ 1.5 mT was applied along the [111] axis 
of the diamond crystal. As explained in the main text, such a magnetic field orientation 
allows both to lift the degeneracy of = ±1 spin sublevels and to divide the NV center 
ensemble into two sub-groups of crystallographic orientations which experience different 
Zeeman splittings. Optically detected ESR spectra were recorded by sweeping the frequency 
of the microwave field while monitoring the PL intensity. As shown in Fig. [5](b), when 
the microwave frequency is resonant with a transition linking = and = ±1 spin 
sublevels, ESR is evidenced as a dip of the PL signal owing to spin-dependent PL response 
of the NV center [1^ 2j. In the following, we focus the study on the ESR line at ujj^i/2ti = 
2.915 GHz, which corresponds to the = ^ = +1 spin transition for the subset 
of [lll]-oriented NV centers. The nitrogen atom of NV centers in our sample being a ^^N 
isotope (99.6% abundance), corresponding to a nuclear spin / = 1, each electron spin state 
is further split into three sublevels by hyperfine interaction with a splitting Ahf = —2.16 
MHz between ESR frequencies associated with different nuclear spin projections [6j. This 
hyperfine structure can be easily observed in our sample by decreasing the microwave power 
in order to reduce power broadening of the ESR linewidth [3j (Fig. [5](c)). 

To probe coherence properties of this subset of NV centers, Ramsey fringes were first 
recorded by using the usual sequence consisting of two microwave tt / 2-pulses separated by a 
variable free evolution duration r (Fig. [5](d)) [HE]. The free induction decay signal exhibits 
beating frequencies which correspond to the hyperfine components of the NV center. Data 
fitting of the free induction decay signal leads to a dephasing time T2 = 390 ± 30 ns of 
the NV center electron spins and its Fourier transform spectrum reveals the ^^N hyperfine 
structure with a linewidth (FWHM) F = 810 ±90 kHz for each peak (see inset of Fig.[5](d)). 

The dephasing time can be greatly enhanced by decoupling the electron spin from its local 
environment with a spin echo sequence (Fig. [5](e)). Using this technique, the dephasing time 
of the NV center ensemble reaches T2 = 7.3 ± 0.4 /xs at room temperature. At high spin 
densities, this quantity is limited by the interaction with a bath of paramagnetic impurities 
including NV centers themselves and PI centers [7^ 8j. 
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B. Superconducting circuit fabrication and measurement setup 

The superconducting circuit is fabricated on a silicon chip oxidized over 50 nm. A 150 nm 
thick niobium layer is first deposited by magnetron sputtering and then dry-etched in a 
SFq plasma to pattern the readout resonator i?, the bus resonator 5, the current lines for 
frequency tuning, and the input waveguides. Finally, the transmon qubit Q, the coupling 
capacitance between Q and S, the Josephson junction of i?, the SQUID in 5, are fabricated 
by double-angle evaporation of aluminum through a shadow mask patterned using e-beam 
lithography. The first layer of aluminum is oxidized in a — O2 mixture to form the oxide 
barrier of the junctions. The chip is glued with wax on a printed circuit board (PCB) and 
wire bonded to it. The PCB is then screwed in a copper box anchored to the cold plate of a 
cryogen-free dilution refrigerator. A complete scheme of the measurement setup and fridge 
wiring is shown in Fig. [6j 

C. Qubit readout 

The qubit readout method we use is explained in detail in [U]. It relies on the nonlinearity 
of the readout resonator R operated in the so-called JBA mode where it behaves as a sample- 
and-hold detector. More precisely, we apply a readout pulse of frequency (jj/27t = 3.456 GHz 
slightly lower than the resonance frequency (jJr/27t = 3.468 GHz, and of power Pr chosen 
so that the resonator is driven close to its bistability, in a regime where the field inside 
the resonator can switch from a low- amplitude state L to a high- amplitude state H. This 
switching is easily detected by measuring the phase of the refiected readout pulse. Repeating 
the same sequence then yields the resonator switching probability Pgw for a given readout 
pulse power. This allows to reconstruct so-called S-curves Psw{Pr) which change from 
to 1 in a narrow power range, close to bistability (see Fig. [7|. Due to the qubit-resonator 
dispersive coupling, the readout resonator R frequency is shifted by a qubit-state-dependent 
quantity ±x so that for well-chosen pulse frequency and power the qubit state is mapped 
onto the resonator dynamical state at the end of the readout pulse (see Fig. [7|. The 
measured switching probability Pgw is therefore directly linked to the qubit excited state 
probability P^ which is the quantity of interest in our experiment. 

That Pg is not directly given by Pgw is due to readout errors caused either by a too small 
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Figure 6: Measurement setup and wiring, a, Full configuration of the measurement apparatus 
at room temperature, b, Scheme of the wiring inside the dilution refrigerator. LPFl, LPF2, 
LPF3, and LPF4 are low-pass filters with cutoff frequencies 1.35 GHz, 450 MHz, 5.4 GHz, and 5.4 
GHz, respectively. BPFl is a band-pass filter with a bandwidth of 2.5 - 4 GHz. CuNi coax is a 
coaxial cable made of CuNi, and Ag-CuNi coax is a silver-plated CuNi coaxial cable. SC coax is 
a superconducting NbTi coaxial cable. Flex coax is a low-loss flexible coaxial cable. Rectangles 
represent ports terminated by 50 ft. The cryogenic microwave amplifier is a CITCRYO 1-12 from 
Caltech, with gain ^ 38 dB and noise temperature ^ 5 K at 3 GHz. A DC magnetic field Bj^y is 
applied parallel to the chip by passing a DC current through an outer superconducting coil. The 
sample box and the coil are surrounded by two magnetic shieldings consisting of a lead cylinder 
around which permalloy tape is wrapped. The sample box, coil, and the shieldings are thermally 
anchored at the mixing chamber. 



X or by qubit relaxation between the end of the experimental sequence and the time at 
which readout effectively takes place. These errors can be modelled with two parameters : 
the probability Cq that the resonator switches despite the qubit being in \g) at the end of 
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Figure 7: Readout resonator switching probability Pgw as a function of readout pulse power 
(S-curves). Black solid line : S-curve with the qubit in thermal equilibrium. Olive solid line: S- 
curve with the qubit prepared in |e) by a tt pulse just before the readout. Red solid line: S-curve 
with the qubit prepared in |e) by a tt pulse, measured with a composite readout pulse including a 
TT pulse on the |e) — |/) transition followed by the usual readout pulse. Dashed lines represent fits 
of these S-curves using a sum of three Erf functions (corresponding to the three transmon states 
1^) , |e) , I/)) with different weights. In this way the thermal population of the qubit Pe,eq — 0.08 
is evaluated as explained in the text. Note that this thermal excitation is responsible for the "foot" 
of the black curve. Dotted blue vertical line indicates the readout power used for simple readout 
pulses, dotted brown vertical line indicates the readout power used for the composite readout pulse 
method. 

the experimental sequence, and the probability ei that the resonator doesn't switch while 
the qubit is in |e). In order to determine eg and ei, we measure the switching probability 
Pswo for a qubit at thermal equilibrium, and after a tt pulse Pswtt that we assume ideal in 
the sense that it swaps states l^') and |e) with 100% efficiency. An additional complication 
arises from the fact that the qubit has a small but finite probability Pe^eq to be found 
in |e) even at thermal equilibrium, due to the rather low qubit frequency chosen in the 
experiment to match the NV centers. We therefore first estimate P^^eq by fitting the shape 
of S-curves at equilibrium and after a tt pulse to a simple model, yielding P^^eq = 0.08 in our 
experiment (see Fig. [?]). This corresponds to an effective electromagnetic temperature of 
50 mK, slightly higher than the cryostat base temperature 30 mK possibly due to imperfect 
filtering of the flux lines. We then find eg and ei by solving the system of two equations 
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Figure 8: Rabi oscillations. The switching probability Psw of readout resonator R is plotted as a 
function of microwave pulse duration At. Inset: pulse sequence used in this measurement. Each 
data point of Psw was constructed by repeating this sequence 10^ times. 

^5^0 = eo(l - Pe,eq) + (l " ei)Pe,eg and Ps^tt = ^^Pe.eq + (l " ei)(l - Pe.eq) - TMs alloWS tO 

determine Pg from the directly measured Pgw since Pgw = eo(l — Pe) + (1 ~ ei)Pe- 

An additional complication arises from the fact that the fidelity of the readout can be 
enhanced (i.e. eg and ei lowered) by using the second excited state |/) of the transmon: for 
that, one applies a tt pulse on the |e) — |/) transition just prior to readout, resulting in a 
so-called composite readout pulse. As explained in [9j this reduces readout errors caused by 
relaxation during the readout pule. Due to technical complications, we use the composite 
readout pulse method only in experiments reported in figures 2 and 4a of the main article. 
The other experiments were performed with simple readout pulses. As a result two different 
sets of errors Cq and ei were determined for each of the two types of readout pulses. For 
composite readout pulses, we find eo = and eg = 0.1, indicating a very high fidelity readout 
consistent with [9j. Without the composite readout pulse we find Cq = and eg = 0.33. 
From these values we convert the measured Pgw into P^ in all our experiments. 

D. Qubit state manipulation 

Single-qubit operations are carried out by applying Gaussian shaped microwave pulses [TO] 
at (jUq. These pulses are generated as explained above by mixing a CW source at (jUq—uus with 
a Gaussian shaped pulse modulated at us using an IQ mixer. Before the measurement such 
as shown in the main text, the system was calibrated to compensate the mixer imperfections 
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Interaction time, x (ns) Adiabatic ramp duration, x (ns) 

Figure 9: Comparison between resonant and adiabatic SWAP pulses. Left panel: principle of a 
resonant SWAP. After excitation of the qubit in |e), 5 is put suddenly into resonance with Q for 
an interaction time r during which |e, 0) and 1) exchange periodically energy. After a time 
T = 37 ns, the qubit excitation is transferred to B. Right panel: principle of an adiabatic SWAP 
(aSWAP). The qubit is excited in |e), after what cj^ is ramped through ujq in a time r, and the 
state of Q is finally read-out. For long enough ramp durations (for this sequence r > 300 ns), the 
qubit excited state population is fully transferred into the bus. 

(amplitude and phase imbalance, offsets). By changing the sideband frequency ujs it is also 
possible to apply pulses on the |e) — |/) transition as requested sometimes for readout. 
Resulting Rabi oscillations are shown in [8| Here the tt pulse and 7r/2 pulse are defined to 
be 50 ns and 25 ns respectively. 

E. Adiabatic SWAP pulses 

SWAP operations between the qubit and the resonator quantum bus can be performed 
by tuning suddenly cj^ in resonance with ujq for a duration Ti/2gQ. We show in Fig. [ojthe 
resulting vacuum Rabi oscillations. In the experiment however, we found out that such a 
resonant SWAP operation was not stable enough to allow subsequent data acquisition longer 
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than ^ 15 minutes. The problem is caused by flux noise in the SQUID loop which causes 
(jJb{^) to change over time so that the amplitude of flux pulse needed to perform the vacuum 
Rabi oscillations in resonance also changes in time. We note that we found a much larger 
flux noise for B^v = 1-4 mT than B^v = niT, probably due to vortices being trapped in 
the superconducting thin fllms around the SQUID. 

To circumvent this problem, we resort instead to adiabatic SWAP operations in which 
ujb is adiabatically ramped through resonance with ujq so that state |e,0) is adiabatically 
converted into l^', 1) for a sufliciently slow flux pulse, yielding the same operation as the reso- 
nant SWAP (see Fig. |9]). Finding good parameters for the pulse requires some optimization 
since a too fast pulse will not be adiabatic while a too slow pulse will strongly reduce the 
signal because of energy relaxation either in the qubit or in the resonator bus. The flnal 
parameters that we used are: ujb starts at 2.52 GHz, is flrst ramped up to 2.589 GHz in 
60 ns, then to 2.643 GHz in 350 ns, then to 2.687 GHz in 40 ns (with a qubit frequency 
ujQ/27r = 2.607 GHz). 

F. Theory 

We now explain in more detail how the theory curves in flgures 2 and 4 of the article 
are calculated. The calculation assumes that the qubit state is perfectly transferred to the 
resonator bus 5, so that the measured perfectly maps the flnal resonator population in 
the state. Each result of the calculation is rescaled in amplitude and offset to match the 
experimental data (this accounts for the additional losses caused by relaxation of the qubit 
or resonator during the pulse sequence, and in particular during the two aSWAPs). Apart 
from that, all the calculations are performed using the following theory, and one unique 
parameters set for the whole paper. 

1. Rabi oscillations (figure 2) 

What needs to be calculated is the flnal probability to flnd a photon in the bus resonator 
after the resonator-spins interaction, assuming the resonator is in at time r = 0. The 
calculations are performed in the Hoist ein-Primakoff approximation, in which the spins and 
the resonator are described by harmonic oscillators, as explained in the Methods section. 
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The system Hamiltonian is H/h = (jJb{^)(i^cl + Xl^j^j^j + ^^^^(^l^ ~ ^j^^)' 9i being the 
couphng constant of spin j with the resonator. We need to calculate p{t) = \(y{t)f with 
a{t) = (0| a{t)a'^{0) |0), which represents the probability that a photon created at t = is 
still present at time t. As shown in [llj this quantity can be calculated by considering an 
effective non-Hermitian Hamiltonian 



Heff/h 
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with complex angular frequencies Cjb = ujb — if^/"^ and Cok = cjj^ — i7o/2 ; here, 70 is 
the spontaneous emission rate of each spin, and k = ujb/Q the bus resonator damping 
rate (where we introduced its quality factor Q). Indeed, introducing the vector X{t) of 
coordinates [(a(t)a'''(0)) , (6j(t)a'''(0)) , ...] it can be shown that dX/dt = —{i/h)HeffX. 
The formal solution to this equation is then 

X{t)=C-'[{s + iH,ff/h)-'xm, (2) 

with X(0) = xg and xq = (1, 0, 0) . This implies that a{t) = xq^ • X{t) = C'^ ^\{s)\ 
with ^1(5) = xq^ ' (5 + iHeff)~^ • Xg and C[f{s)] = f e~^^ f{t)dt^ s being a complex num- 
ber. Since ^1(5) is not singular on its imaginary axis, we only need ti for pure imaginary 
argument s = —iuj to perfom the Laplace transform inversion. As shown in [llj, we have 
ti{-iuj) = i/ [(jj - (jJb + - W{uj)] with W{uj) = g'^ f p{uj')duj' / [uj - uj' + i7o/2]. In this 
last formula, g is the coupling constant of the spin ensemble to the bus resonator and p{uj) 
is the density of spins which is taken as explained below. Computing a{t) is thus achieved 
by evaluating ti for the distribution p{uj)^ and numerically evaluting the inverse Laplace 
transform. At the end of the calculation, we take the 70 ^ limit since NV centers at low 
temperature have negligible energy relaxation. 

2. Single-photon Ramsey experiment (figure 4b) 

For the Ramsey-like experiment (figure 4b), each 7r/2 pulse is realised by bringing the 
resonator and spins to resonance. For a fast pulse, the resonant interaction maps continu- 
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ously the coherent state of the field to the superradiant mode in the spins. We cahbrate the 
interaction time in such a way as to transform the state xq into the superposition 
After the first 7r/2 pulse, the resonator is kept detuned from the spin ensemble for a time 
t. The system state at this point can be evaluated using eq. (2). We define Xcit) (resp. 
Xs{t)) as the vector of coordinates \(^a{t)a} {Q)) , (hj{t)a^ {Q)) , ...] at time t with initial 
conditions xq (resp. xs)- A second 7r/2 pulse is then applied before the amplitude a{t) of 
the field in the resonator is measured: 

a{t) = l^xl-U^/2{XG{t)-Xs{t)) = ^{xl+xl)iXG{t)-Xs{t)) = ^C-\h{s)-t2{s)+h{s)-t 

(3) 

where t2{s) = xs'^ • {s + iH^ff)'^ • xs and t2{s) = xs'^ • {s + iH^ff)'^ • x^, ti and are defined 
above. As shown in pjj, t2{—iuj) = —ti{—i(ju)W{(ju){s + iuj^)/ and ^3 = — ^4. 

3. Parameters used in the simulation 

We thus see that the only thing that is needed to numerically perform these calculations 
is the density of spins p{uj) and the ensemble coupling constant g. The spin density p{uj) is 
chosen as the sum of three Lorentzian peaks separated by 2.3 MHz. The only parameters 
that are adjusted in order to fit the data are the coupling constants and the peak linewidth 
in the HF structure. As explained in the main text, we find that gijlii = 2.9 MHz and 
5'////27r = 3.8 MHz fit best our data, with a linewidth of 1.6 MHz for the spins belonging 
to group /, and 2.4 MHz for the spins belonging to group ///. We attribute the larger 
linewidth of the spins from group /// to a residual misalignment of B^y with respect to 
the [1,1,1] axis of the crystal which causes each of the three < 1,1,1 > axes non-collinear 
with the field to undergo slightly different Zeeman shifts. A misalignment of 0.02 rad would 
be enough to cause a broadening such as we observe. 

We also note that the splitting of 2.3 MHz between the three peaks of the HF structure 
is slightly larger than the value reported in most articles which is 2.18 MHz. Our data are 
however not precise enough to determine precisely whether this difference actually refiects 
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a change in the HF interaction parameters of the NV center at low temperature. 
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